Stability of 3D Cubic Fixed Point in 
Two-Coupling-Constant ^'^-Theory 

H. Kleinert, S. Thorns and V. Schulte-Frohlinde 
Institut fiir Theoretische Physik 
Freie Universitat Berlin 
Arnimallee 14 
14195 Berlin 

Abstract 

For an anisotropic euclidean (p^-theovy with two interactions 
['^iJ2i=i (Pi)'^ + 4'i] the /3-functions are calculated from five-loop 

perturbation expansions in d = 4 — e dimensions, using the knowledge 
of the large-order behavior and Borel transformations. For e = 1, an 
infrared stable cubic fixed point for M > 3 is found, implying that the 
critical exponents in the magnetic phase transition of real crystals are 
of the cubic universality class. There were previous indications of the 
stability based either on lower-loop expansions or on less reliable Fade 
approximations, but only the evidence presented in this work seems 
to be sufficently convincing to draw this conclusion. 

1 Introduction 



The most elegant approach to phase transitions in many physical systems 
proceeds via field-theoretic renormalization group techniques Jll. The best- 
studied system is the isotropic Heisenberg ferromagnet with M classical spin 
components. Its critical behavior can be described correctly by an 0(M)- 
symmetric vector field theory with a quartic interaction Y^fii 

In a real crystal, such an interaction is never presentably by itself. The 
crystalline structure gives rise to anisotropics, most prominently of cubic 
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symmetry, which can be represented by an extra field interaction Y^fLi (pt- 
This term breaks the 0(M)-symmetry by favoring magnetizations along the 
edges or the diagonals of a hypercube in M dimensions. The extended theory 
interpolates between an 0(M)-symmetric and a cubic system. It has been 
pointed out a long time ago that, depending on M, the 0(M)-symmetric 
and the cubic fixed point interchange their stability. For M < M^, the 
0(M)-symmetric, isotropic fixed point is stable. For M > M^. the isotropy 
is destabilized and the trajectories of renormalization fiow cross over to the 
cubic fixed point. Estimates using calculations up to three-loops in- 
dicated that Mc must lie somewhere between 3 and 4. Resummation proce- 
dures based on Fade approximations P| suggested Mc to lie below 3, thus 
permitting real crystals to exhibit critical exponents of the cubic univer- 
sality class. The uncertainty of these estimates have prompted Kleinert and 
Schulte-Frohlinde |Q to carry the expansions up to five loops. They increased 
the evidence for M^. < 3 considerably, again via Fade resummation. 

For a simple (^^-theory, the Fade approximation is known to be inaccurate. 
At present, the most accurate renormalization group functions for that theory 
have been obtained by combining perturbation expansions with large-order 
estimates, using a resummation procedure based on Borel-transformations 
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Intending the application of these more powerful resummation methods, 
the large-order behavior of renormalization group functions has recently been 
derived for the M- vector model with cubic anisotropy, by Kleinert and Thoms 



It is the purpose of this letter to combine these large-order results with 
the five-loop perturbation expansions [^] using a simple Borel-type of resum- 
mation algorithm 0, whose power has been exhibited in recent model 



studies by Kleinert, Thoms and Janke |]T5|. 

The results to be presented in this paper allow us to conclude with a 
reasonable certainty that an infrared-stable cubic fixed point exists at the 
physically most relevant value M = 3. However, due to the vicinity of the 
isotropic fixed point, the differences in the critical exponents are very hard 
to measure experimentally. Going beyond the Fade work in Refs. [|] and 
we also show explicitly the instability and stability of the isotropic and the 
cubic fixed point, respectively. 
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2 Resummation 



2.1 The Problem 

Object of investigation is a 0^-theory with cubic anisotropy. The correspond- 
ing energy functional reads: 



1 Svr f Ub ^ Vb 



^df,(j)Bid^(j)Bi + — iy — Sijkl + -^^ijkl ) (l)Bi(j)Bj(f)Bk(j)Bl 



where (pBii^) {i = 1,2,..., M) is the bare M-component field in d = A — e 
dimensions, and ub, vb are the bare coupling constants. In particular, we 
shall consider the physically most interesting case of M = 3, and continue e 
to e = 1. The tensors associated with the two interaction terms in (|l|) have 
the following symmetrized form: 

Sijkl = -{^ij^kl + ^ik^jl + ^il^jk) , 



5, 



i = j = k = I 
'^''^ ~ 1 , otherwise . 



The symmetry of the action under reflection 0j —>■ —(pi and under permuta- 
tions of the M field indices i implies the following form of the vertex functions, 
persisting to all orders in perturbation theory: 

rSf ~ ; vf^ ~ r(^'^)5., , (2) 

^%-^^S,,ki + T^^^,H. (3) 

The symmetry permits us to renormalize the M components of the field 
and the composite field with only a single renormalization constant 
and Z02, respectively. The bare field 0b, the composite field |0^, and the 
two coupling constants ub and vb are related to the corresponding physical 
objects by 



(x) = iZ,, 



/^2 



(X), 



Ub = ^fZuiZ^y^u] Vb = ^fZ^{Z^)-'^v 

where /i is a mass parameter. We employ dimensional regularization with 
minimal subtraction. The square brackets around 0^ indicate a renormaliza- 
tion of this composite operator. Recall that this renormalization is different 
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from that of the wave function which converts 0b to 0. In fact it is closely 
related to the mass renormalization of a theory with mass m 7^ 0. In the 
minimal subtraction scheme, the corresponding renormalization constants 
are related by: Z^2 = Z„2(Z0)~^. 

The RG-functions are defined in the usual way: 



(5) 



The natural parameter for the anisotropy of the system is the ratio 6 = 
v/{u + v), and the isotropic case corresponds to 5 = 0. We shall use the new 
couplings g = u + V and 6 for the calculation of the fixed points from the 
resummed /3-functions 

(3<^{g, 6) = f^'^Hg, 6),v{g, 6)] + f3^[u{g, 6),v{g, 6)] , 
gP\g,6) = {1 ~ 6)(3^[u{g,6),v{g,6)]- 6p-[u{g,6),v{g,6)]. (6) 

The 0(M) -symmetric and cubic fixed points can be obtained by calculating 
the simultaneous zeros {g*,6*) of P^[u{g,S),v{g,6)] and (3^[u{g,6),v{g,6)]. 
For the physically interesting number of field components, M = 3, the 
infrared- stable cubic fixed point is expected to appear very close to the 0(M)- 
symmetric one. Since 6 is very small in this region, it will be sufficient to 
restrict the resummation efforts to the (7-series accompanying each power 5", 
so that the /3-functions at the cubic fixed point will be approximated by 

N 

= nuig\6*),vig\6*)] ^ J2 K^''\9*)S*\ 

n=0 

N 

= P^Hg*, n, v{g\ 5*)] ^ 5* B<''\g*)5*''-' , (7) 

n=l 

where B^"* (g) = res Y.k=n Pkng^ indicates resummed g-seiies. 

From the five-loop perturbation expansion in Ref. [^], the perturbation 
coefficients f3^l^^ are known up to the order N = 6. For e = 1 and the number 
of field components M = 3, the following expansions are known: 

6 

^ P^^gf' = ~g + 3M7g^ - 7M7g^ + 47.651^^ - 437.646^^ + 4998.62c/^ , 

A:=0 
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^ f31^^g>' = g- 5.333C/2 + 15.667^^ - 121.767/ + 1341.05^^ - 17821.1/ , 

k=l 
6 

Pl^gk = 1.667^2 - 10/ + 115.885/ - 1664.86c/^ + 27191/ , 

k=2 
6 

Pl^g^ = 2/ - 50.074/ + 1064.62c/^ - 22916.2/ , 

k=3 
6 

Y Pl^g^ = 8.305/ - 350.528^^ + 11183.1/ , 

fc=4 

Y iild = 47.368^7^ - 2966.14^6 , (3^,g' = 330.76^^ , 

k=5 

and 

6 

Y Pkid'' = -9 + 4/ - 10.778^^ + 75.875/ - 776.26c/^ + 9707.36^^ , 

k=l 
6 

Y 131^9^ = -g^ + 6.222/ - 67.319/ + 944.05/ - 15030.9c/^ , 

k=2 
6 

J2 Pk39'' = -1-111/ + 30.211/ - 639.243^^ + 13549.6/ , 

A:=3 
6 

Y I3l^g^ = -6.218/ + 233.262c/^ - 7122Mg^ , 

k=A 

^/?i;5/ = -33.414/ + 1973.58/, /^^V = -228.19/ . (8) 

fc=5 



In addition to the five-loop expansions, the large-order behavior of and 



has been obtained explicitly in Ref. ||1J], with the result 



Pl'n i-l f k\ A;('^+5)/2+" [1 + 0{l/k)] , A; » n . (9) 

The resummation algorithm to be employed in this paper will make use of 
all these informations. 



2.2 The Algorithm 

As explained in detail in Ref. [0, it is possible to reexpand a divergent 
perturbation series 

oo k 

Z{g,5) = Y.Y.Zkng'S-. (10) 

fc=On=0 

in a special infinite set of Borel summable functions Ipn{g) as 



so that the approximation 



n=0 



.p=n 



(11) 



N 



Z{g,6)^Z'^''\g,6) = Y: 



n=0 



N 



^ ] ^pn -l^pn{g) 



p=n 



N 



= Zi''\g)6- (12) 



n=0 



has the same series as Z{g, 6) up to the powers g'^S'^ , while reproducing the 
known large-order behavior of the perturbation expansion ([T0|). 



Zkn ^ l{n) {-ay k\ k^- [1 + 0{l/k)] , k^n. 



(13) 



As shown by Janke and Kleinert in Refs. 0, H, the natural choice for the 
functions Ipn{g) are certain confluent hypergeometric functions, where for the 



case of two coupling constants a second index was introduced in Ref. 
One possibility was to use 



Ipnia) 



(14) 



6o(n)+l 



dw 



(1 + w;)w^o(")+P 



with 



r[6o(n) + 1] (1 - w;)2''o(n)+2a 



&o(^) =/?« + -• 



+3 



exp 



Aw 



\ — w^ag 



(15) 



Then the coefficients apn are given by 



"pn 



E 



■'kn 



n (Mn) + l)fc 



p + k-l-2a 
p — k 



(16) 
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where Ck = r(c + k)/r{c) are Pochhammer's symbols. The parameter a is 
free to choose, and may be used to accommodate any strong-couphng power 
oiZig,S) 

Zig,5)'^ (17) 

if this is known. Since in quantum field theory, this is not the case, a will be 
chosen by the condition of best convergence of the resummed ^f-series Z^^\g) 
in Eq. (]T2|), as explained in the next section. 



2.3 Optimal Choice of Strong-Coupling Power a 

An idea of the relevance of a is gained from the study of corresponding models 
in quantum mechanics, where the strong-coupling behavior can be deduced 
from scaling arguments. Consider first an 0(M) -symmetric anharmonic os- 
cillator with g{x1Y /A interaction, where M is the number of components of 
the vector x. Here the functions zl^\g) in Eq. (|12]) represent the resummed 
ground-state energies E^^\g). In Figure |I] we have illustrated the conver- 
gence of E'^^'^ for the anharmonic oscillator with one x-component (M = 1) 
at a coupling constant g/A = 0.1. We have plotted E^^^ versus the order 
of approximation for various values of the strong-coupling parameter a. 
At large A^, the curves become increasingly independent of a and approach 
a saturation value which coincides with the ground-state energy. This ex- 
act result does not depend on a. Therefore, we choose the strong-coupling 
parameter under the condition that the curvature and the slope of the cor- 
responding curve depend minimally on the variation of a, when approaching 
the saturation region. This choice complies with the principle of minimal 
sensitivity (PMS) which has been used with great success to optimize vari- 
ational perturbation expansions [jl6|. Using the discretized form of the first 
and second derivative, the optimal a-value, to be denoted by OpMS; is found 
by calculating the extrema or turning points of the judicial function 

A(^=)(^7,a) = (18) 

where A"s indicates the beginning of the saturation region. In our special 
example, = 7 (see Figure ^a). In Eq. (|TBp the value of the coupling 
constant g is chosen such that the error of the resummation becomes small 
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and the influence of a is isolated. Since the rate of convergence of the re- 
summation decreases for increasing we determine the optimal a at a small 
coupling constant gjA^ = 0.1. In Figure |T].6, we have plotted the a-dependence 
of A'^=((yf, a) for Ng = 7 and g/4 = 0.1. The optimal a-value is found to be 
•^PMS ~ 0.341, which is very close to the exactly known strong-coupling pa- 
rameter a = 1/3. 

We mention that for the zero-dimensional case which corresponds to a 
simple integral with 5fx^/4-interaction, the above criterion yields the exact 
strong-coupling parameter a = —1/4 (see Figure |^). 

In most applications of quantum field theory, perturbation series are too 
short to detect the formation of the saturation plateau with sufficient accu- 
racy. We shall see in the following that for shorter series the criterion given 
above can be simplified without a significant loss of accuracy by neglecting 
the curvature term in Eq. (|TBp. Then the slope term is deduced from the 
two highest known partial sums i^;(^max-l) g^j^^ ^(A^max) -^y forming: 

where it is important that _£;{^max-i) g^^^ ^(A^max) j-^ave gotten over the 
large-fiuctuating initial region. Thus the optimal a-value, denoted now by 
«PMS; will be determined from calculating the extrema or turning points of 
^(Afmax)('(^^ Q,^y'^(Afmax-i)^^^ Q,')^ beforc, thc coupllug constant is chosen to 
be small (7/4 = 0.1 to ensure a sensitive determination of the optimal a. 

As the largest available order of approximation for the RG-functions in 
QFT is iVmax = 6, we have plotted in Figure ^ the a-dependence of the 
ratio E^^^ / E^^^ for the numbers of vector components M = 1 and M = 2 
at a coupling constant g/A = 0.1. In the case of M = 1, a minimum exists 
at apMS = 0.389. For M = 2, there is no extremum, and the optimal a- 
value lies at the turning point apMS = 0.323. In both cases the simplified 
criterion yields again results for apMS which are close to the exactly known 
strong-coupling parameter a = 1/3. 

It is unnecessary to know a to a higher accuracy than that. In order to 
show this, we have compared the resummed ground-state energies E(^) of 
the anharmonic oscillator using various values of a (see Table |1|). For the 
coupling constants g/i = 0.1 and g/A = 1.0, the rate of convergence to the 
exact energies depends only very little on a. 
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Consider now the case that Z{g, 6) in Eq. (|1^) represents the ground-state 
energy E{g,6) of the anisotropic quartic oscillator with an interaction 



x^ + 2{l-S)x^y^ + y^ (20) 



mt ^ 

which we have studied in detail in [|l^ using the Borel-type resummation 
algorithm of Refs. §]. When expressed in terms of the old coupling 
constants u, v via g = u+v and 6 = v/ (u+v), the expression (^) corresponds 
to the interaction term in the euclidean action (|1|) for the number of field 



components M = 2. It has been found in that the parameter a is the 
same for each coefficient En{g) in the (5-expansion of the ground-state energy 
E{g,S) = J2'^^=o ^n{g)S"', and having at each n the same value a = 1/3 
as in the isotropic case. The approximation E^\g) which follows from a 
resummation of the corresponding perturbation series in g up to the order 
is becoming less accurate for increasing n since the gf-expansions have 
fewer and fewer terms. However, taking into account the smallness of the 
anisotropy 5, the 5-expansion may be truncated at a finite order A^. This 
yields the approximation 

N 
n=0 

which was found to be very accurate in a wide region of 5 around 5 = 0. We 
have compared the result for the ground-state energy resummed at a = 1/3, 
which is known from Ref. [|l^], with the result obtained for apMS = 0.323. 
Figure || shows the ^-dependence of the approximated ground-state energy 
E^^\g, 5) for the two different values of the parameter a and various coupling 
constants g/A. For g/A = 0.1, the two curves for a and apMS coincide. From 
Figure § and Table |I] we can thus conclude that the error which is caused by 
an inaccurate determination of a is negligible. 

2.4 Application to Quantum Field Theory 

The above analysis of the anisotropic oscillator can now be applied to the 
corresponding model in quantum field theory. The resummation of the (3- 
functions is carried out by combining the formulas ([T^) , (|14D, and ([T6|) , 



where the function Z stands now for /3'^ and P"". The parameters bo{n) and 
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a follow from the large-order behavior and are the same for both (3- 
functions: 

3 

ho{n) = I3n + - = Q + n, 

a = 1. (21) 

The optimal value of the parameter a is chosen to cause minimal sensitivity 
of the ratios B^^^^ / Bf^ and B^^/Bf^ on a at the small value of coupling 
constant g/A = 0.1 [recall Eq. (0)] in accordance with the above observations 
for the 0(M)-symmetric anharmonic oscillator. 
The optimal values apMs are found to be: 

apMs = 1.348 , P'' : apMs = 1-225 . (22) 

For the simultaneous solution of Eqs. (|^), we have first determined for each 
/9-function all zero-point functions 5^''^ (g) which are implicitly defined by 

/?"/l^?,^iX(^)] = 0. (23) 

From the second equation in (^, we have then read off a trivial solution 
5^^\g) = 0. Restricting attention to the region around the isotropic limit 
5 = 0, we have found numerically a second nontrivial solution S^^\g). For /?" 
only one solution 5„((?) was found. Having obtained 6!^^^ and 6u the isotropic 
fixed point follows from the condition 

'^:iso(^?-o) = Clsi(^?L) = 0. (24) 

The cubic fixed point is similarly obtained by calculating the solution of the 
equation: 

^ucub(5'cub) = '^?)cub(fi'cub) ■ (25) 

3 Results 

Table ^ contains the numerical values of the isotropic and the cubic fixed 
points for increasing order of approximation A^. Starting from the order 
= 3 a cubic fixed point is obtained which lies in the upper half of the 
plane of coupling constants u and v {6 > 0). Plotting the functions Su{g) 



10 



and 5lj^\g), yields the cubic fixed point (fi'cub) '^cub) the crossing point. 
This is shown in Figure ^ for the orders of approximation N = 2 and N = 6. 
The cubic fixed point is found to he very close to the isotropic one. 

In order to convince ourselves of the stability of the cubic fixed point at 
the number of field components M = 3, we calculate the eigenvalues bi and 
62 of the matrix 

B = ( r:. "z. 1 (26) 




using the resummed /^-functions (|]). The result is contained in Table ^ 
If both eigenvalues are positive, the corresponding fixed point is infrared 
stable. For M = 3, this is definitely the case for the cubic fixed point. At the 
isotropic fixed point, one the other hand, one eigenvalue 62 is negative. As 
the isotropic and the cubic fixed point interchange their stability at M = Mc, 
the result corroborate the suggestion in ||^ and [H that the critical value Mc 
lies below M = 3. 

Thus we conclude that the critical behavior of magnetic phase transitions 
in anisotropic crystals with cubic symmetry is governed by the cubic, not by 
the isotropic Heisenberg fixed point. The corresponding critical exponents r] 
and V follow from the resummed RG-functions 7^ and 7,^2 via the defining 
relations 

r/ = 70(^7*,5*), u-'-2 = ~^^.{g\5*). (27) 

Unfortunately, our result is only of fundamental interest and has no easily 
measurable experimental consequences. Due to the vicinity of the isotropic 
fixed point, the difference in the critical exponents is smaller than one per- 
cent, so that the new cubic universality class is practically indistinguishable 
from the isotropic class. 
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Figure Captions 



Figure 1: Convergence of the ground-state energy E of the anharmonic os- 
cillator with (yf/4 = 0.1. a) Resummed ground-state energies are plotted 
versus the order of approximation N for various strong-coupling parameters 
a. h) The a-dependence of the function A*^^) of Eq. (p^S]). The optimal 
a-value is given by the minimum ttp^^g = 0.3408, very close to the exactly 
known value 1/3. 

Figure 2: Judicial function A'^''^ of Eq. (|TB|) to determine the strong- 
coupling power a for the simple integral Z{g) = J dxexp{—x'^ /2 — gx'^/4), 
plotted as a function of a for g/A = 0.1. The optimal parameter a^j^ lies 
at the lowest value of A^'^\g,a), and is found to be equal to the exact one. 
It can be shown that for a = Oexact = — 1/4, already the zeroth order of the 
resummation algorithm used in this paper reproduces the exact value of Z{g) 
This is the origin of the cusp and A'^^''\g, a) = at a = —1/4. 

Figure 3: The a-dependence of the ratio E'-^^ / E^^'^ in Eq. (||) for the 0(M)- 
symmetric anharmonic oscillator at constant coupling strength g/A = 0.1 
and various numbers of field components M. a) For the simple anharmonic 
oscillator (M = 1) the optimal a- value apMS = 0.389 lies at the minimum, b) 
For the two-component oscillator (M = 2) the optimal value apMS = 0.323 
lies at the turning point. 

Figure 4: Ground state energy E of the anisotropic oscillator with the 
interacting potential (PP]), as a function of the anisotropy parameter 6 for 
two coupling constants g/4 = 0.1 and g/A = 1.0. Comparison is made 
between the approximation E^^^ obtained once for the exact strong-coupling 
parameter Oexact = 1/3 and once for apMS = 0.323. In the case g/i = 
0.1, differences are invisible at this graphical resolution. For comparison, 
results from another resummation scheme, variational perturbation theory 
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(VPT) obtained in Ref. |T^ are shown as well. 



Figure 5: Determination of the cubic fixed point g* , 5* . The nontrivial 
zeros 5u and of the /3-functions /3" and (3"" are plotted against g [see Eqs. 
(0) and (p3|) ]. The the values of g* and 5* are found from the intersection 
point, a) Up to the order of approximation = 2, no cubic fixed point exist. 
h) For the order of approximation = 6, the cubic fixed point is given by 
g* = 0.39154, 6* = 0.015309. 
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Figure 1: Convergence of the ground-state energy E of the anharmonic os- 
cillator with g/A = 0.1. a) Resummed ground-state energies are plotted 
versus the order of approximation N for various strong-coupling parameters 
a. h) The a-dependence of the function A*^^) of Eq. (p^S]). The optimal 
a-value is given by the minimum apMs = 0.3408, very close to the exactly 
known value 1/3. 
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Figure 2: Judicial function /S.'^'^\g^a) of Eq. (|18]) to determine the strong- 
coupling power a for the simple integral Z{g) = J (ia:exp(— — gx'^/A), 
plotted as a function of a for g/A = 0.1. The optimal parameter a^jjj lies 
at the lowest value of A^'^\g,a), and is found to be equal to the exact one. 
It can be shown that for a = aexact = —1/4, already the zeroth order of the 
resummation algorithm used in this paper reproduces the exact value of Z[g) 
[^]. This is the origin of the cusp and A^'^''\g, a) = at a = —1/4. 
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Figure 3: The a-dependence of the ratio E^^^ E^^'^ in Eq. (|19|) for the 0(M)- 
symmetric anharmonic oscillator at constant coupling strength g/A = 0.1 
and various numbers of field components M. a) For the simple anharmonic 
oscillator (M = 1) the optimal a- value apMs = 0.389 lies at the minimum, b) 
For the two-component oscillator (M = 2) the optimal value apMs = 0.323 
lies at the turning point. 
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Figure 4: Ground state energy E of the anisotropic oscillator with the 
interacting potential (PP]), as a function of the anisotropy parameter 5 for 
two coupling constants g/A = 0.1 and g/A = 1.0. Comparison is made 
between the approximation E^^^ obtained once for the exact strong-coupling 
parameter Oexact = 1/3 and once for apMs = 0.323. In the case g/A = 
0.1, differences are invisible at this graphical resolution. For comparison, 
results from another resummation scheme, variational perturbation theory 
(VPT) obtained in Ref. [T^ are shown as well. 
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Figure 5: Determination of the cubic fixed point g* ^ 6*. The nontrivial 
zeros 5u and 5':jp of the /3-functions (3^ and (3^ are plotted against g [see Eqs. 
(|^ and (P^l)]. The the values of g* and 5* are found from the intersection 
point, a) Up to the order of approximation = 2, no cubic fixed point exist. 
h) For the order of approximation = 6, the cubic fixed point is given by 
g* = 0.39154, 6* = 0.015309. 
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77T 


g/4 = 0.1 


g/4 = 1.0 


N 


Q^exact OipMS 


Q^exact OipMS 


1 

2 

3 
4 
5 
6 
7 
8 
9 


0.561496 0.56235 
0.558592 0.558614 

0.559232 0.559254 
0.559142 0.559142 
0.559143 0.559143 
0.559147 0.559147 
0.559146 0.559146 
0.559146 0.559146 
0.559146 0.559146 


0.83055 0.849631 
0.78297 0.784942 

0.812948 0.816638 
0.801761 0.802012 
0.802206 0.802487 
0.805103 0.805518 
0.803901 0.803937 
0.803115 0.803072 
0.803852 0.803924 


exact 


0.559146 


0.803770 



Table 1: Convergence of the ground-state energy Eq of the anharmonic 
oscillator with M = 1 for the strong-coupling parameters ccexact = 1/3 and 
QipMS — 0.389. 



N 


q* 

JISO 


5.* 

ISO 


fl'cub 


"cub 


2 


0.560616 





don't exist 


3 


0.440796 





0.50208 


0.291074 


4 


0.393506 





0.400199 


0.037862 


5 


0.4012 





0.411057 


0.063068 


6 


0.389037 





0.39154 


0.015309 



Table 2: Numerical result for the isotropic and the cubic fixed point for 
increasing order of approximation N. For iV > 3, a cubic fixed point is 
found in the upper half of the coupling constant plane {u, v), i. e. 6*^^ > 0. 
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l,cub 




LISO 


4 


0.782796 


0.0048920 


0.784532 


-0.00502046 


5 


0.764835 


0.00851725 


0.763966 


-0.00886277 


6 


0.80609 


0.00212717 


0.80658 


-0.00214788 



Table 3: Stability of cubic fixed point, as demonstrated by the eigenvalues of 
the stability matrix B in Eq. (|26|) , calculated from the resummed /3-functions 
Pg and Ps in Eq- (|§)- The isotropic fixed point is unstable. 
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